Non- asymptotic fractional order differentiators via 
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Abstract — Recently, Mboup, Join and Fliess |27|, |28| intro- 
duced non-asymptotic integer order differentiators by using an 
algebraic parametric estimation method |7|, |8|. In this paper, in 
order to obtain non-asymptotic fractional order differentiators 
we apply this algebraic parametric method to truncated expan- 
sions of fractional Taylor series based on the Jumarie's modified 
Riemann-Liouville derivative 1 14 1. Exact and simple formulae for 
these differentiators are given where a sliding integration window 
of a noisy signal involving Jacobi polynomials is used without 
complex mathematical deduction. The efficiency and the stability 
with respect to corrupting noises of the proposed fractional order 
differentiators are shown in numerical simulations. 

I. INTRODUCTION 

Fractional models arise in many practical situations (J3), 
(371 for example). Such fractional order systems may also be 
used for control purposes: CRONE control is known to have 
good robustness properties (see l33ll . (34], ll35l . 11371 ). In order 
to implement such controller one needs to have a good digital 
fractional order differentiator from noisy signals, which is the 
scope of this paper. 

The fractional derivative has a long history and has often ap- 
peared in science, engineering and finance (see, e.g., Il32l . [36|, 
IE], 0). Different from classical integer order derivative, 
there are several kinds of definitions for the fractional deriva- 
tive which are generally not equivalent with each other |36|, 
[ IT | . Among these definitions, the Riemann-Liouville deriva- 
tive and the Caputo derivative are often used [36 1. Recently, 
a new modified Riemann-Liouville derivative is proposed by 
Jumarie B141 . This new definition of fractional derivative 
has two main advantages: firstly comparing with the Caputo 
derivative, the function to be differentiated is not necessarily 
differentiable, secondly different from the Riemann-Liouville 
derivative the Jumarie's modified Riemann-Liouville derivative 
of a constant is defined to zero. Moreover, a fractional Taylor 
series expansion |[T3l was invented by using this new defini- 
tion. Thanks to these merits, the Jumarie's modified Riemann- 
Liouville derivative was successfully applied (see, e.g., [15|, 
ifTBI . |03]). The fractional order differentiator is concerned 
with estimating the fractional order derivatives of an unknown 
signal from its noisy observed data. Because of its importance, 
various methods have been developed during the last years 



(ED, [25 1 for example). However, the obtained fractional order 
differentiators were usually based on the Riemann-Liouville 
derivative and the Caputo derivative. To our knowledge, there 
is no one based on the Jumarie's modified Riemann-Liouville 
derivative. 

Recent algebraic parametric estimation method for linear 
systems Q, 0, (44ll has been extended to various problems 
in signal processing (see, e.g., Il9l. iflOl. 11261. l30l. iTJTI. |42|. 
[41 1, [19|, [18]). Let us emphasize that this method is algebraic 
and non-asymptotic, which provides explicit formulae and 
finite-time estimates. Moreover, it exhibits good robustness 
properties with respect to corrupting noises, without the need 
of knowing their statistical properties (see 0, for more 
theoretical details). The robustness properties have already 
been confirmed by numerous computer simulations and several 
laboratory experiments. Very recently, this method was used to 
solve the ill-posed numerical differentiation problem in (27], 
[28 1 . The main idea is to apply a differentiator called integral - 
annihilator to a truncated Taylor series expansion which is a 
local approximation of the signal to be differentiated. Stable 
non-asymptotic differentiators of integer order were exactly 
given by the integrals of noisy signals involving the Jacobi 
orthogonal polynomials. The associated estimation errors were 
studied in [20|, [22], ||231 . An extension of these differentiators 
for multivariate numerical differentiation were proposed in 
[ 38 1, 1 39 1 , l40l . However, this method has not been used to 
estimate fractional order derivatives. 

The aim of this paper is to introduce non-asymptotic 
fractional order differentiators for the Jumarie's modified 
Riemann-Liouville derivative by using the algebraic paramet- 
ric method. In Section [II] we recall how to apply the algebraic 
parametric method to obtain an integer order differentiator 
from a truncated Taylor series expansion. Section [TIT] be- 
gins with the definition of the Jumarie's modified Riemann- 
Liouville derivative. Then, two fractional order differentiators 
are obtained by applying integral-annihilators to truncated 
fractional Taylor series expansions with different truncated 
orders. Moreover, it is shown that the differentiator obtained 
from higher order truncated expansion can be expressed as 
an affine combination of the ones obtained from lower order 



truncated expansion. Numerical tests are given in Section |IV] 
They help us to show the efficiency and the stability of the 
proposed fractional order differentiators. Finally, we give some 
conclusions and perspectives for our future work in SectionlVl 

II. METHODOLOGY 

Let y = x + xn be a noisy signal observed in an open 
interval Jcl, where x G C (I) with neN and the nois^H 
w is bounded and integrable. In this section, we are going to 
recall how to use an algebraic parametric method to estimate 
the integer order derivatives of x [27|, |28|. 

For any to g /, we introduce the set D t(l := {f £ R+; to + 
t G /}. By using the famous Taylor's formula formulated by 
Hardy (HI] p. 293), we obtain that Vt 6 J, t € D to , 



" fj 



ast^O. (1) 



In a similar way to classical numerical differentiation 
methods, the algebraic parametric method uses the n th order 
derivative of a polynomial to estimate the one of x. Here, we 
consider the following truncated Taylor series expansion of x 
on M+ 



Vt G 



x n (to + 1) 



E^ (j) ('o)- 



(2) 



Different from the other methods, the n th order derivative of 
the approximation polynomial x n is calculated by applying 
algebraic manipulations to x n in the operational domain. 
Precisely, we apply a differential operator with the following 
form El) 



n 



1 



(3) 



where — 1 < /i G M and /j G N. This differential operator was 
introduced in |28] with /i G N. 

Let us mention that the operation £ n+k ■ s" is used to 
annihilate all the terms containing x^(to) with < j < n— 1 
in x„ in the operational domain, and the rational term 
permits to obtain a Riemann-Liouville integral ll24l expression 
for x^{to) in the time domain (see the next section for more 
details). Thus, the differential operator nj! is called integral- 
annihilator for x^ n \to) via x n ESI . Finally, if we replace x n 
by y in the obtained integral expression, then we can obtain the 
following estimator for x^(to) (see [22 1 and [28 1 for more 
details) 



Vt G/, y t ( n) (fc,/i,T) := 

^ln,k,» J q w^ k {r)P^ k \r) y(to + Tr) dr, 



(4) 



'More generally, the noise is a stochastic process, which is bounded with 
certain probability and integrable in the sense of convergence in mean square 
(see (SI). 



where T G D u> , = B(n+fc+i,n+ M +i) ' B (''') is the 

classical beta function ((TJ, p. 258), Pn 1 is the n th order 
Jacobi polynomial ([1 1 p. 775) defined on [0, 1] as follows 



j=0 



n + [i, 
j 



k 



n- 3 



: Ut-1T- j (t)\ (5) 



and w^ t k{T~) = (1 — r) AI r is the associated weight function. 
Hence, this differentiator depends on three parameters fc, /i, T. 
Moreover, it is a non-asymptotic pointwise differentiator using 
the sliding integration window [to,to + T]. Since y[™\k, n,T) 
is obtained by taking an n th order polynomial where n is the 
order of the derivative estimated, we call it minimal Jacobi 
differentiator. 

III. NON-ASYMPTOTIC FRACTIONAL ORDER 
DIFFERENTIATOR 

In this section, we are going to propose two fractional order 
differentiators by applying the algebraic parametric method to 
different truncated expansions of the fractional Taylor series 
based on the Jumarie's modified Riemann-Liouville derivative 

m, ma. 

A. The Jumarie's modified Riemann-Liouville derivative 

Let / be a continuous function defined on M, then the Ju- 
marie's modified Riemann-Liouville derivative of / is defined 
as follows lfl4l 

1 d' r< f{r) - /(0) 



/W(t) := 



T(l-a)dt l Jo (t-r) a+1 



■ dr, 



(6) 



where < I — 1 < a < I with I e N*. This fractional order 
derivative is in fact defined through the fractional difference 



(7) 



where h > 0, FW f(t) = fit + h) and 
A a f(t) = (FW-l) a f(t) 

=£(-!)*(?)/ [*+(«-<)*]. (8) 

i=o v J 

Let us recall that the Riemann-Liouville derivative is defined 
as follows ((36) p. 62) 



m 

T{l-a)dt l J ( t - T ) 



where < I - 1 < a < I with / G N*. If we take f(t) = f ' 
with n G N and t G MI, then we obtain (see O p. 72) 



T(n + 1) 
r(n + l-a) 



with a G 



(10) 



Consequently, the Jumarie's modified Riemann-Liouville 
derivative can be expressed as follows 



Vt G R+, f (a \t) = D?f(t) 



r 



T(l-a) 



/(o). 



(11) 



One of some useful properties of the Jumarie's modified 
Riemann-Liouville derivative is the fractional Leibniz deriva- 
tive rule lfT4l 



(fg) {a) = f {a) g + g (a) f 



(12) 



Moreover, since x G C n (I), a generalized Taylor series 
expansion based on the Jumarie's modified Riemann-Liouville 
derivative of x is given as follows [14|, [13]: \/t G /, t G D to , 



f^nn + n + i) 



x^ +n \t ), 



(13) 

where 7 = a — n with n < a < n + 1. 

Let us remark that the fractional order derivative x^ 7+ ™- ) 
in ( [T3] l should be understood as (x^)^' 1 ^ which is different 
from (x^)( n > lfT4l . 

B. Minimal Jacobi fractional differentiator 

By using ([T3~l >, we take the following truncated fractional 
Taylor series expansion of x on M + : Vto £ I, t £ M + , 



c Q (to + t) := ^ ^a:«) (t ) + w ^ ^ (*o ) ■ 'I^S 
j=o J ' 



r(a+l) 



Then, by using the algebraic parametric method we can give 
the following proposition. 

Proposition 1: Let y = x+zu be a noisy signal observed on 
an open interval Jcl, where x G C n (I) with n < a < n+l, 
n G N, and w be a bounded and integrable noise. Then an 
estimator for the a th order derivative value x^ a '(to) is given 
by: Vt G /, 

^(fc^T) := 

:^7n,fc,M,« y w H,k { T ) P n+l ( T ) 2/ (*o + Tt) dr, 



'J'o 



where T G A , 7n , fc , M , a = B(a+ ?£ n "j M+2) , and is 
the Jacobi polynomial defined by <j3j with k G N and — 1 < 

/iel. 

Proof. By applying the Laplace transform to (fT4|) . we get 

n 

i Q ( a ) =^ s -^ +1 )xW(to) + s- (Q+1) xM(^o), (16) 
3=0 

where x a (s) is the Laplace transform of x a (to + 1), and s is 
the Laplace variable. 

We are going to apply some algebraic manipulations to ( fToT l. 
Firstly, we apply the operation ^„ +1+fc • s n+1 to (IT6b so as 
to annihilate the terms containing x^-'(io) with < j < n in 
x Q . Thus, we get 

jn+l+fc / 2 



=( _ ir+1+ «i>±i± j e_j_ l( „, ((0) . 



Secondly, if we apply the Leibniz derivative rule to 
x a (s), then the highest order of s in the obtained 



n+l 



ds" + 1 + kS 

sum is n+l. Hence, we choose the rational term S „ + Vn with 



fi > — 1 such that n + 2 +v ■ s 



n+l _ 



x-ftt with fi + 1 > 0. This 



permits to obtain a Riemann-Liouville integral B24II expression 
of x a in the time domain. Consequently, we construct an 
integral-annihilator for x^ a '(to) via x a which corresponds to 
the operator n^ 1 defined by OJ. 

Consequently, by applying the inverse Laplace transform 
and the classical rules of operational calculus, we get 

^ j 11 '^ 1 (j2s-^x^(t ) + s -( Q+1 V")(t )j j (T) = 

( ir+1+fc £(a + l + fc) 2"+°+'+*+" g(a 
r(a — n) r(n + a + k + [i + 3) 

(T - T) n+ ^ +1 T n+k+1 X^ +1) {t + T)dT. 

(17) 



T(n + fjt + 2) 7 



By applying a change of variable r — > TV and n + l times 
integrations by parts, we get 



— ar(«)(to) = 

in,k,/j,,a 

1 

__ {(1 _ r) n+P+l T ™+fe+l} ^ + Tt) dr> 



(18) 



Finally, this proof can be completed by substituting x Q in (fT~8T > 
by y and applying the Rodrigues formula (JT| p. 785) to the 
right side of (T8]l. □ 
Since the differentiator y["\k,fx,T) involves a Jacobi 
polynomial, and it is obtained by taking the minimal order 
(j = 1) truncated fractional Taylor series expansion in ( fT~3b , 
we call it minimal Jacobi fractional differentiator. The corre- 
sponding truncated error part comes from the truncated term 

— — rZ~TT\ x (*<>)• If we take a = n + 1 in (B), 



J=2 



r(j 7 + n + 1) 



then it is easy to obtain that 

ylf(k,n,T) = y[: +1 \k,n,T), 



(19) 



where y^? (k, fi, T) is the minimal Jacobi differentiator 
for x( n+ V(t ) given in @. Moreover, in a similar way to 
the minimal Jacobi differentiator method, the minimal Jacobi 
fractional differentiator y^'(k,fJ,,T) can also be obtained 
by using the classical orthogonal properties of the Jacobi 
polynomials to ( fL4T > (see [21] for more details). Consequently, 
as done in 11221 . the parameter k which is defined on N can 
also be extended to ] — 1, +00 [. 

Finally, let us recall that the affine Jacobi differentiator 
being an affine combination of the minimal Jacobi differen- 
tiators was introduced in ||27l , (28 1 by applying the algebraic 
parametric method to a higher order truncated Taylor series 
expansion than (fJJ. Hence, the convergence rate for the Jacobi 



differentiator was improved. Similarly to the integer derivative X2a-n m the time domain. Consequently, n£"t a is an integral- 
case, we are going to introduce an Affine Jacobi fractional annihilator for x (a) (i ) via x 2a -n- Moreover, we have 

(_!)»+!+* £-1 {ll n k ^ a x 2a _ n (s)} (T) = 
r(a + 1 + k) (a - n) T n + a+3+k +» {a) < 25) 

~n7 \ "FT i i 7 i i \ 0): 



differentiator in the next subsection. 



C. Affine Jacobi fractional differentiator 



In this subsection, we take a new truncated fractional Taylor 
series expansion in £□) up to j = 2: V t £ I, t £ R+ r(a - ra) T(n + a + fc + M + 4) 



,T 2 . 



tj where T £ D to . □ 

q-wfto + 1) = -rrff^fto) Similarly to Proposition Q] by substituting a^a-n 

3=0 (20) by w in the obtained Riemann-Liouville integral 



t 2a " x [2a-n), t \ ^ 1 |n^ Q i 2Q -n} we obtain a new estimator 



r ( a + *) r ( 2a - n + !) for 3fW(i ). We denote it by y t ( a) (fc, /ti, T, a). Then, we get 

In order to obtain an estimator for x^ (t ) from this truncated the following proposition. 

expansion, we introduce the following differential operator Proposition 3: Let y be a noisy signal defined as in Propo- 
sition [T| then we have the following affine relation 

1 tt rjn+k+l 
yjn+1 _ 1 il 2Q-n+l + fe _ „n+l 

s2Q+fe+ , I+3 ds ds n + k + i > y t ( o a) (fc, Ai ,T,a) = A Q , fe ,„y t ( a) (fc, / , + l,T) (26) 

+ (l-A Q , fe , n )y t ( a) (fc + l,M,T), 

where -K/ieR, fceN and n < a < n + 1 with n £ N. 

Then, we give the following proposition. where A Qifei „ = 2a ~"J" ra 1+fc . The differentiator y t ( " } (fc, jii, T, a) 

Proposition 2: The differential operator II£+* Q defined in i s thus called affine Jacobi fractional differentiator. 

d2T} is an integral-annihilator for x( Q) (t ) via a; 2Q! -n defined In order to prove the above proposition, we need the 

(|20). following lemma. 



in 

Proof. Applying the Laplace transform to ([20]), we get Lemma 1: Let / be the Laplace transform of an analytic 

„ function f defined on I, and II? +1 be the integral- 

X2a-n( s ) = ^ s~^ +1 ^ x^\to) annihilators defined in (fJTJ and (01 respectively, then we have 

+ s -(a+ l ) a .(a)/ to) + fl -(2«-n+X) a .(2a-n) (t0 J n £, M ,a/( s ) ~ 

/ .j t \ T-rll 4-1 n / \ -w—w-'il _1_ 1 n / \ 

where X2a-n(s) is the Laplace transform of X2a-n(fo + t). 



(2a-n+l + k) U^ + J(s) + IT^/( S ). 



We are going to apply lH£t a to (1221 . Firstly, we apply the Proof. By applying the Leibniz derivative rule, we get 
operation ^ s ^+i +k ■ s n+1 to (l22l so as to annihilate the terms n" +1 
containing a;^'(to) with < j < n in ^2 a _ n . Thus, we obtain 



1 / Jn+fc+l 
1 4.*.r(a+ 1 + fc) 1 ,„ w . + - 2a-n+k+ia n+1 f (\ 



ds n+fc+2 ' 
rn+1 , Tl n + 1 



n+1+k r(2a-n + l + k)x^(t ) =(2a-n + l + fc)n^ +1 /( S )+n^/( S ). 
+( r( 2a _ 2n) s 2 Q -™+fc+i ■ 

Secondly, we apply the operator £ ■ s 2 Q -«+i+fe t0 gjj so Proof of Proposition H According to <ES>, we obtain 

as to annihilate the term containing x^ 2a ^ n " > (to) in a^a-n- _i f „ + i , s~l 

This yields " £ \ n fc,i*,a £ 2a-n(s)) (T) 



□ 



" + ' +/ - S " +1 (i2a-n( S ))= =£_1 C E , " (j ' +lll(<1 ( i «) + , " t " +1,l(0l ( 1 »)| 



+1+fc r(a + 1 + fc) _ n)sa -»- lx(a)(to)- 

r ("-«) = £ (nfc^^^i 5 )) ( T ) 



Thirdly, if we apply the Leibniz derivative rule to 



_ ( ir+1+fe r(a + 1 + k) (a - n) T"+"+3+fe+M ^ 



d_ g 2 a "+ 1 + fc ^ +1+fc S "+ 1 x 2a ,»(s), then the highest order ; ' ; T(a-n) T(n + a + k + p, + 4)' 

of s in the obtained sum is 2a + n + k + 2. Hence, we (28) 

apply the rational term sststetstt; with /i > — 1 such that 

l _ g 2a+n+k+2 _ ^^^^ _|_ > q -p n j s By substituting x a by a; in the Riemann-Liouville integrals 

allows us to obtain a Riemann-Liouville integral expression of £■ i^jfe^a^aC 8 )} (-0 ^ 1 i^feji ^a( s )| (^) gi ven i n 



and ( fT7T > respectively, we get 

4 Q) (fc,M,T,a) : =5h^l £ -i {n^f (a)} (T) 
4 Q) (fc, /x + 1,T) :=c niklP , a) r£- 1 {n'fc++i^(^)} (r) 

+ 1 )M ,T) := - ^f^-£ 1 jllSi^W} CO, 



(29) 



whprp r , ( l~\n+l+k r(a-n) r(«+a+fc+p+4) 

wnere c n ^,fi, a ,T — { J-) r(a+i+*0 T"+°+ 3 + fc +" ■ 

Hence, by using Lemma [T] we obtain 



^ \k,n,T,a) = A Q , fe ,„xJ" ; (fc,^ + 1,T) 

+ (l-A Q , fe ,„)4 Q) (fc + l^,7 1 )- 



(30) 



Finally, this proof can be completed by substituting x by y in 

4 a) (fc, M) T, a), (fc + 1, /i, T, a) and z£> (fc, /x + 1, T, a). 

□ 



Noisy observation and original signal 




Fig. 1. Signal x and noisy signal y. 




2 3 
Time 



Fig. 2. y t (0 - 5) (0, 0, T) with T = 0.25 and y t (0 5) (0, 0, T, 0.5) with T - 
0.26. 



IV. NUMERICAL SIMULATIONS 

In order to show the efficiency and the stability of the pro- 
posed fractional order differentiators, we give some numerical 
results in this section. 




Fig. 3. yf- 7) (0, 0, T) with T = 0.25 and ^°' 7) (0, 0, T, 0.7) with T 
0.28. 



We assume that j/(tj) = + is the discrete 

noisy observation of x on I = [0,4] where x(ti) = 
exp(0.2tj) sin(5tj), zu is a zero-mean white Gaussian noise, 
and U = T s i for i = 0, ■ ■ ■ ,4 X 10 3 with T s = ^ The 
variance of is adjusted in such a way that the signal- 
to-noise ratio SNR = 101og 10 ( ^^(top ) is e q ual t0 
SNR — 28.07dB. We can see the original signal x and 
its noisy observation y in Figure [1] The exact Jumarie's 
modified Riemann-Liouville derivative of x can be calculated 
by using ( fT2l . ( fTTT i and the Riemann-Liouville derivative of 
the functions sin(-) and exp(-) (see ||29l , p. 83). We estimate 



the derivatives with a — 0.5 and 0.7 respectively. 

We apply the trapezoidal numerical integration method to 
approximate the integrals in our differentiators where we use 
m + 1 (T = mTs) discrete observation values in each sliding 
integration window. 

The formal derivatives and their estimated values are shown 
in Figure |2] and Figure [3] where the parameters k and ^ are set 
to zero for the Jacobi fractional differentiators. On one hand, 
the error due to the noise for the minimal Jacobi fractional 
differentiator can be negligible with respect to the truncated 
error part which produces a time-shift in the estimate. On the 
other hand, the noise error for the affine Jacobi fractional 
differentiator is larger than the one for the minimal Jacobi 
fractional differentiator, but the truncated error part is much 
smaller. Hence, the time-shift is significant reduced. 

V. CONCLUSION 

In this paper, two non-asymptotic fractional order differen- 
tiators called minimal and affine Jacobi fractional differentia- 
tors are proposed by applying an algebraic parametric method 
to truncated expansions of fractional Taylor series. They can 
be used to estimate the recently invented Jumarie's modified 
Riemann-Liouville derivative. Numerical simulations are given 
so as to show their efficiency and stability with respect to 
corrupting noises. It is shown that these differentiators can also 
be obtained by using the classical orthogonal properties of the 
Jacobi polynomials [21 1. By using this method, a generalized 



affine Jacobi fractional differentiator will be obtained from a 
truncated fractional Taylor series expansion with an arbitrary 
truncated order in the future work. Moreover, in a similar way 
to ll2D and ll22l . the noise errors and the truncated errors in 
these differentiators will be analyzed. In particular, we will 
study the influence of the parameters k, /1 and T to these errors 
so as to give a guideline for choosing the optimal parameters. 

References 

[1] M. Abramowitz and LA. Stegun, editeurs. Handbook of mathematical 
functions. GPO, 1965. 

[2] D.L. Chen, Y.Q. Chen and D.Y. Xue: Digital Fractional Order Savitzky- 
Golay Differentiator, IEEE Transactions on Circuits and Systems II: 
Express Briefs, 58, 11, pp. 758-762, 2011. 

[3] M. Cugnet, J. Sabatier, S. Laruelle, S. Grugeon, B. Sahut, A. Oustaloup 
and J.M. Tarascon: On Lead-acid battery resistance and cranking capa- 
bility estimation, IEEE Transactions on Industrial Electronics, 57, 3, pp. 
909-917, 2010. 

[4] M. Fliess and R. Hotzel: Sur les systemes lineaires a derivation non 
entiere. C.R. Acad. Sci. Paris Ser. lib, (signal, informatique), 324, pp. 
99-105, 1997. 

[5] M. Fliess: Analyse non standard du bruit, C.R. Acad. Sci. Paris Ser. I, 
342, pp. 797-802, 2006. 

[6] M. Fliess: Critique du rapport signal a bruit en communications 
numeriques - Questioning the signal to noise ratio in digital commu- 
nications, International Conference in Honor of Claude Lobry, Revue 
africaine d'informatique et de Mathematiques appliquees, 9, pp. 419- 
429, 2008. 

[7] M. Fliess and H. Sira-Ramfrez: An algebraic framework for linear 
identification, ESAIM Control Optim. Calc. Variat., 9, pp. 151-168, 
2003. 

[8] M. Fliess and H. Sira-Ramfrez: Closed-loop parametric identification 
for continuous-time linear systems via new algebraic techniques, in H. 
Gamier, L. Wang (Eds): Identification of Continuous-time Models from 
Sampled Data, pp. 363-391, Springer, 2008. 

[9] M. Fliess, M. Mboup, H. Mounier and H. Sira-Ramfrez: Questioning 
some paradigms of signal processing via concrete examples, in Algebraic 
Methods in Flatness, Signal Processing and State Estimation, H. Sira- 
Rarmrez, G. Silva-Navarro (Eds.), Editiorial Lagares, Mexico, pp. 1-21, 
2003. 

[10] M. Fliess, C. Join, M. Mboup and H. Sira-Ramfrez: Compression 
differentielle de transitoires bruites. C.R. Acad. Sci., 1(339), pp. 821- 
826, Paris, 2004. 

[11] GH. Hardy: A course of pure mathematics. Cambridge University Press, 

tenth edition, 1952. 
[12] R. Hotzel and M. Fliess: On linear systems with a fractional derivation: 

Introductory theory with examples. Mathematics and computers in 

Simulation, 45, pp. 385-395, 1998. 
[13] G. Jumarie: Modified Riemann-Liouville derivative and fractional Taylor 

series of nondifferentiable functions further results, Comput. Math. Appl. 

51, pp. 1367-1376, 2006. 
[14] G. Jumarie: Table of some basic fractional calculus formulae derived 

from a modified Riemann-Liouville derivative for non-differentiable 

functions, Appl. Math. Lett. 22, pp. 378-385 2009. 
[15] G. Jumarie: Fourier's transform of fractional order via Mittag-Leffler 

function and modified Riemann-Liouville derivative, J. Appl. Math. 

Inform. 26, pp. 1101-1121, 2008. 
[16] G. Jumarie: Laplac's transform of fractional order via the Mittag-Leffler 

function and modified Riemann-Liouville derivative, Appl. Math. Lett. 

22, pp. 1659-1664, 2009. 
[17] A. A. Kilbas, H.M. Srivastava, and J.J. Trujillo: Theory and Applica- 
tions of Fractional Differential Equations, vol. 204 of North-Holland 

Mathematics Studies, Elsevier, Amsterdam, The Netherlands, 2006. 
[18] D.Y. Liu, O. Gibaru, W. Perruquetti, M. Fliess and M. Mboup: An error 

analysis in the algebraic estimation of a noisy sinusoidal signal. In: 

16th Mediterranean conference on Control and automation (MED'08), 

Ajaccio, France, 2008. 
[19] D.Y. Liu, O. Gibaru and W. Perruquetti: Parameters estimation of a noisy 

sinusoidal signal with time-varying amplitude. In: 19th Mediterranean 

conference on Control and automation (MED' 11), Corfu, Greece, 2011. 



[20] D.Y. Liu, O. Gibaru and W. Perruquetti: Error analysis for a class 
of numerical differentiator: application to state observation, 48th IEEE 
Conference on Decision and Control, Shanghai, China, 2009. 

[21] D.Y. Liu, O. Gibaru and W. Perruquetti: Differentiation by integration 
with Jacobi polynomials. J. Comput. Appl. Math., 235, 9, pp. 3015-3032, 
2011. 

[22] D.Y. Liu, O. Gibaru and W. Perruquetti: Error analysis of Jacobi 
derivative estimators for noisy signals. Numerical Algorithms, 58, 1, 
pp. 53-83, 2011. 

[23] D.Y. Liu, O. Gibaru and W. Perruquetti: Convergence Rate of the Causal 

Jacobi Derivative Estimator. Curves and Surfaces 2011, LNCS 6920 

proceedings, pp. 45-55, 2011. 
[24] A. Loverro: Fractional calculus, history, definitions and applications for 

the engineer. Rapport technique, Univeristy of Notre Dame: Department 

of Aerospace and Mechanical Engineering, May 2004. 
[25] J.A.T. Machado: Calculation of fractional derivatives of noisy data with 

genetic algorithms, Nonlinear Dyn., 57, 1/2, pp. 253-260, Jul. 2009. 
[26] M. Mboup: Parameter estimation for signals described by differential 

equations, Applicable Analysis, 88, pp. 29-52, 2009. 
[27] M. Mboup, C. Join and M. Fliess: A revised look at numerical 

differentiation with an application to nonlinear feedback control. In: 

15th Mediterranean conference on Control and automation (MED'07). 

Athenes, Greece, 2007. 
[28] M. Mboup, C. Join and M. Fliess: Numerical differentiation with 

annihilators in noisy environment, Numerical Algorithms 50, 4, pp. 439- 

467, 2009. 

[29] K.S. Miller and B. Ross, An Introduction to the Fractional Calculus and 
Fractional Differential Equations, Wiley, New York, 1933. 

[30] A. Neves, M. Mboup and M. Fliess: An Algebraic Receiver for Full 
Response CPM Demodulation, VI International Telecommunications 
Symposium (ITS2006), September 3-6, Fortaleza-ce, Brazil, 2006. 

[31] A. Neves, M.D. Miranda and M. Mboup: Algebraic parameter estimation 
of damped exponentials, Proc. 15 th Europ. Signal Processing Conf. - 
EUSIPCO 2007, Poznari, 2007. 

[32] K.B. Oldham and J. Spanier: The Fractional Calculus, Academic Press, 
New York 1974. 

[33] A. Oustaloup, B. Mathieu and P. Lanusse: The CRONE control of 

resonant plants : application to a flexible transmission, European Journal 

of Control, 1, 2, pp. 113-121, 1995. 
[34] A. Oustaloup, J. Sabatier and X. Moreau: From fractal robustness to the 

CRONE approach, ESAIM: Proceedings, pp. 177-192, Dec. 1998. 
[35] A. Oustaloup, J. Sabatier and P. Lanusse: From fractal robustness to the 

CRONE control, FCAA, 1, 2, pp. 1-30, Jan. 1999. 
[36] I. Podlubny: Fractional Differential Equations, vol. 198 of Mathematics 

in Science and Engineering, Academic Press, New York, NY, USA, 

1999. 

[37] V. Pommier, J. Sabatier, P. Lanusse and A. Oustaloup: CRONE control 
of a nonlinear hydraulic actuator, Control Engineering Practice, 10, pp. 
391-402, Jan. 2002. 

[38] S. Riachy, M. Mboup and J. P. Richard: Multivariate numerical differen- 
tiation. J. Comput. Appl. Math., 236, 6, pp. 1069-1089, 2011. 

[39] S. Riachy, Y. Bachalany, M. Mboup, and J. P. Richard: An algebraic 
method for multi-dimensional derivative estimation. In: 16th Mediter- 
ranean conference on Control and automation (MED'08), Ajaccio, 
France, 2008. 

[40] S. Riachy, Y. Bachalany, M. Mboup and J. P. Richard: Differenciation 
numerique multivariable I: estimateurs algebriques et structure, in: 
6 ,eme Confrence Internationale Francophone d'Automatique, 2010. 

[41] J.R. Trapero, H. Sira-Ramfrez and V.F. Battle: An algebraic frequency 
estimator for a biased and noisy sinusoidal signal, Signal Processing, 
87, pp. 1188-1201, 2007. 

[42] R. Ushirobira, W. Perruquetti, M. Mboup and M. Fliess: Estimation 
algebrique des parametres intrinseques d'un signal sinusoidal biaise en 
environnement bruite, XXIII e Colloque Gretsi, Bordeaux, France, 2011. 

[43] GC. Wu and E.W.M. Lee: Fractional Variational Iteration method and 
Its Appliation, Phys. Lett. A 374, pp. 2506- 2509, 2010. 

[44] Y. Tian, T. Floquet and W. Perruquetti, Fast state estimation in linear 
time-varying systems: an algebraic approach, 47th IEEE Conference on 
Decision and Control, Cancun, Mexique, 2008. 



Estimations 




